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Abstract 

Spatial models for spread of an epidemic may be mapped onto bond percolation. 
We point out that with disorder in the strength of contacts between individuals 
patchiness in the spread of the epidemic is very likely, and the criterion for epi- 
demic outbreak depends strongly on the disorder because the critical region of the 
corresponding percolation model is broadened. In some networks the percolation 
threshold is zero if another kind of disorder is present, namely divergent fluctua- 
tions in the number of contacts. We give an example, a network with a well defined 
geography, where this is not necessarily so, and discuss whether real infection net- 
works are likely to have this property. 
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1 Introduction 

Models for the spread of epidemics are an interesting application of non- 
equilibrium statistical physics. In this paper we discuss two aspects of this 
sort of study in cases when disorder is present. First we generalize the well- 
known [1] application of percolation theory to epidemics for which variation 
in the strength of coupling. Second, we consider the contact networks which 
carry epidemics. If there is large variation in the number of contacts between 
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members of the population, then the network structure can lead to a percola- 
tion threshold of zero, for example if the distribution of contacts is a suitable 
power law. We show an example of a network which seems to have Pc 7^ for 
the same power-law distribution, and which might have some relevance to real 
infections. 



2 SIR epidemics in an inhomogeneous population 

The SIR (Susceptible, Infected, Recovered) model is a classic way to think 

about the spead of epidemic infections [3]. It is based on the idea that in 
many epidemics an individual catches a disease, infects others at rate x for 
a certain period, r, and then recovers or dies, and is thus removed from the 
infection network. A mean-field version of this model is given by the SIR 
equations[3] for a perfectly mixed population: 



Here N is the size of the population, and we normalize x in this way so that we 
can think about the infection rate per contact. Clearly, we have an outbreak 



We can think of Eq. (1) as describing dynamics on a fully conjugated graph - 
every node sees every other equally. In this paper we will consider a model on 
a lattice, where each individual only sees its neighbors. Each site of the lattice 
is initially occupied by a susceptible (S) individual. One lattice site is initially 
infected, /. Any S can be infected by Ps which are nearest neighbors. The 
probability of infection per unit time along the i^^ bond, Xi, is chosen from a 
distribution f{x). For a case near the threshold for the epidemic (see below) 
the pattern of infection is shown in Figure 1. The similarity of this pattern to 
a percolation cluster is very clear. 

The existence of a large epidemic depends on the parameters just as in the 
mean-field case. For fixed Xi, if the recovery time, r, is too small the epidemic 
will die out, and if r is large enough it will persist. A critical value, Tc marks 
the threshold. However, even for large r there is a finite probability that the 
epidemic dies. We define a spanning probability, M(r), the probability that 
an infection started from a single site in the middle of a large lattice reaches 
the edge. M grows rapidly for r > Tc- In Figure 2 we show M for a 200 x 200 
lattice as a function of r for two different choices of f{x), one rectangular, and 
one strongly asymmetric. 

This model is not obviously related to percolation theory. For a very narrow 
distribution, f{x) = S{x — Xg) Grassberger [1] showed that we do have perco- 



dS/dt 



xSI/N, 



dl/dt = xSI/N -I/t 




if Ro = xS{0)t/N > 1. 
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Fig. 1. Pattern of recovered after the 
epidemic has died out on a 256x256 
lattice. Recovered sites are gray, and 
uninfected black. 




Fig. 2. M(r) for a triangular 
lattice, A, and a square lattice, 
□ . Solid lines:/(x) from Eq. (3) 
Xmax = 0.003. Dotted lines: Eq. (4) 

with Xm.in. - 
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lation as follows: the probability to infect a neighbor is p{t) = 1 — (1 — XoY- 
Thus p is the fraction of bonds completed in a given epidemic, and the infec- 
tion process is mapped onto bond percolation with p(r) playing the role of 
the percolation probability. For a square lattice in two dimensions, Pc = 1/2, 
so we expect that when p is greater than 1/2 we will have an epidemic. This 
gives a critical recovery time: Tc = — ln2/ln(l — Xo)- For any other network, 
Pc = 1 — (1 — XoY" determines the threshold for epidemic spread. For arbitrary 
/'s, we can write: 



p{t) 



1 - 



xY f{x)dx = Pc 



(2) 



Thus, for any / we have mapped our infection problem to percolation. The 
solution of Eq. (2) gives Tc. The mapping of a dynamic problem onto static 
percolation is not trivial: different runs for infection spreading in a system with 
a given realization of Xi give different, but not independent realizations of a 
percolation problem, i.e. of a set of a completed bonds. In previous work [5] 
we have verified Eq. (2) by doing a data collapse for various different choices 
of f{x). After data collapse, the function M, expressed as a function of p using 
Eq. (2) is identical to the mass of the infinite cluster in ordinary percolation 
theory. All of the results of this work may be applied easily for any network. 

In order to understand the effect of the distribution of bond strengths we 
studied /'s which are either weakly or strongly asymmetric. An example of a 
symmetric / is: 



fw (■^) 1/ -^n 



< X < Xr 



(3) 
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In cases like this where f{x) is concentrated near its mean value, it is easy to 
show the spread of Xi does influence the infection propagation very much. 

However, the dependence p(t), Eq.(2), is strongly nonlinear and is dominated 
by the behavior of f{x) in vicinity of x — 0: the spread of the infection 
is controlled by weakest bonds. Thus strongly asymmetric distributions con- 
centrated near zero should be very different from the case discussed above. 
Consider, for example, power-laws f{x) = {a — l)x^"' with < a < 1. In 
this case /(I — xY f{x)dx = (a — 1)-B(t -|- 1, 1 — a), where B{x,y) is the 
beta-function. Strongly asymmetric distributions correspond to values of a in 
vicinity of 1. Using the fact that for such distributions r ^ 1, we can use the 
Stirling formula to get in leading order: Tc — [r(2 — q;)(1 — Pc)]~^^''^~"^ which 
diverges for a — > 1. 

An even more extreme case is that of functions f{x) with a > 1. These are 
not acceptable probability densities since they are not normalizable on [0,1]. 
However, one can truncate the / in the vicinity of zero to allow normalization. 
For example: 

fs{x)^C/x Xmin<X<l; C = 1/| ln(x„j„) | . (4) 

The two distributions considered in Figure 2 are given by Eq. (4) with Xmin — 
e~^^, an extremely asymmetric case, and Eq.(3), the symmetric case. 

In Figure 2 we see that the effect of strong asymmetry is to spread out the 
critical region, i.e. the region where M is not close to either 1 or 0. This is 
the parameter region where the epidemic pattern is patchy, as in Figure 1. 
For a strongly asymmetric / the patchy pattern corresponds to a large range 
of the parameter r. We recall from percolation theory that the regime where 
M is neither very small nor very close to 1 is where the fractal nature of the 
percolation cluster at Pc persists over large length scales. This could have im- 
plications for practical epidemiology: if the actual distribution of Xi is very 
asymmetric, we would be very likely to observe patchiness. For such distri- 
butions mean-field theory is very inaccurate [5] and sampling of populations 
must be done with care because of the spatial correlations. 



3 The network of contacts: power-law distributions 

In the previous section we considered networks for which the number of con- 
tacts between individuals does not vary much, though wc supposed that the 
strength of such contacts varied a good deal. However, there is data to sug- 
gest [6] that some infection networks have a wide distribution of the number 
of contacts, so that the bond number. A;, is drawn from a distribution of the 
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form P{k) oc with 2 < a < 3. In this case the mean value of the number 
of contacts, < k >, exists, but the fluctuation, < > diverges. 

We have seen that it is sufficient to consider percolation on such networks 
to understand SIR epidemics. This problem has been discussed a good deal 
[7,8,9]. There is general agreement that pc is zero for many cases. The graphs 
considered by these authors are essentially random graphs with a given degree 
distribution. Thus the model has a kind of random mixing. 

The random graphs dominate the the ensemble of all graphs with a given 
P{k). The study of this ensemble was introduced by Molloy and Reed [2] and 
a criterion for the formation of a giant component was derived. Callaway, et. 
al [8] and Cohen, et al. [7] express a related result for percolation: Pc(< > 
/ < A; >^ —2) = 1 Thus if < A;^ > diverges we expect pc — 0. We can see 
this by starting at a random node on the graph so that there are ko bonds of 
which p/co are active. If neighbors are chosen at random, the probability for 
the number of bonds on the next neighbor is p'^kiP{ki) / < k >. The factor 
of ki comes from the fact that highly bonded sites are more likely to have 
a bond with the first site. Now the number of infected increases by a factor 
< ki > / < ko >— p < k^ > / < k >^ at each such step. If < A;^ > is finite, 
then this can be made less than unity for p small, and the infection will die 
out. U < k^ > is infinite, then p must be 0. 

Now consider a model which has the same P{k) but has geography [10,11]. 
We assume that individuals infect their neighbors, but that the number of 
infections is drawn from P. We construct this by taking each point on a d- 
dimensional lattice with probability p and connecting them all other points 
within a d-sphere of volume k, centered on the point. In 2 dimensions we sur- 
round each point with a disk of radius of roughly R = [k/irY^'^. The probability 
distribution of R is given by P{R) = P{k)dk/dR oc R-^da-d+i) _ Qur model 
is a variant of an old two-dimensional circle model of continuum percolation 
[12] with variable radius of circles. Sec Figure 3. This is a small part of the 
MoUoy-Reed ensemble, but we suggest that it may be relevant for real-world 
epidemics: sometimes infection is primarily local. 

Note that a site can have the bonds it generates itself ('proper' bonds) and 
those from other disks. We show that the total number of bonds per node 
has the same power-law behavior as the proper bonds if a > 2 by calculat- 
ing the mean number of bonds connecting a node i to nodes j whose cen- 
ter is a distance r > Ri away. The probability to find Rj larger than r is 
/^°° P{r')dr'. Thus, the mean number of nodes connected to i from outside is 
kout drr'^~^ P{Rj)dRj. Thus kout — Ri "^^ which tends to zero for 

large R, as long as a > 2. This means that the probability distribution of the 
number of the 'proper' bonds and the actual number of the bonds of a node 
is the same for k large, and that bonds from outside are rare. 
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Fig. 3. The geographic modeh One disk with radius Ra is centered around node A, 
and another with radius i?^ is centered about B. All nodes within each disk are 
attached to the central node. 

For (i = 1 it is easy to see that the system never percolates, because for any 
p < 1 there will eventually be a gap of length A in the chain of connected 
disks. Consider disks with radii between k and k + dk; the concentration of 
their centers is pP{k)dk. The probability that the gap is not covered by one of 
these disks is equal to the probability that none of their centers are found in 
the interval of length 2k + A (in units of the lattice constant) centered at the 
gap, 1 — pP{k){2k + /S)dk ~ exp[— pP(A;)(2A; + /S)dk] The overall probability 
that no disk overlaps the gap is the product over k: 

exp ^- JpP{k) {2k + A)dk^ = exp {-2p (k) - pA) , (5) 



which is finite as long as (k) exists. 

In 2 dimensions we did simulations to find pc- Disks were randomly added to 
an L X L lattice at different sites until a cluster spanned the lattice. Figure 4 
shows the average percolation threshold Pc{L) as a function of It appears 
from the plot that there is a finite percolation threshold for not just a > 3 but 
2 < a < 3 as well, the region of interest for real world epidemics. For a < 2, 
the results are consistent with pc = and seem to scale according as L^"^""^) 
for sufficiently large L. We have shown [10] that this scaling is related to the 
presence of giant disks that span the whole sample. These occur with finite 
probability for a < 2. 

We can try to understand the results in Figure 4 by calculating < ki > as 
above. Most of the bonds emanating from the center of disk point within it. 
These have degrees drawn from the distribution P{ki), not kiP{ki). Thus, 
for the typical behavior we < ki >=< k > for a > 2. However, there are 
rare events which correspond to the overlap of large disks from outside with 
the center of the starting disk. These give a large contribution. In fact, in 2 
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Fig. 4. pc on an L X L lattice as a function of 1/L. From top to bottom, the data 
are for a = 2.5, 2.3, 2.05 (dotted lines), 2, (solid line), and 1.95, 1.9, and 1.7 (dashed 
lines). The solid straight lines have slope 2{a — 2). 

dimensions for a disk of size Rq we add up the contributions of outside disks 
of radius Ri. 

oo oo oo 

< A;i >oc J rdr J RlP{Ri)dRi = J drr'^^''-^^ (6) 

Ro r Ro 



This diverges for a < 3. Note that this is an overestimate since have not 
attempted to exclude already infected sites. 

We have a confusing situation: the typical behavior is that the cluster of 
disks is finite in extent, but the average mass added diverges. There is a 
rigorous theorem in the study of continuum percolation which treats exactly 
this situation [12]. It states that if < k'^~^/'^ > is finite, but < k"^ > diverges 
( in 2 dimensions 2.5 < a < 3), then there is a Pc such that for p < 
with probability unity all clusters are finite, but the expected number of points 
of the largest cluster is infinite. (The probability to span a large lattice is 
not addressed by this theorem.) We interpret this odd situation, at least for 
a > 2.5 as a form of intermittancy. Typical epidemics are finite for small p, but 
huge, rare events occur. It is interesting to ask the physical question of what 
we would expect to observe in this case. We claim (perhaps too hopefully) 
that the situation is best represented by the computer simulations. Figure 4 
shows that at a = 2.5 there is quite solid evidence for a non-zero pc, and this 
seems to extend down to a = 2. (Note that the theorem does not exclude this 
possibility). 
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In summary we have a counterexample to the prevaihng notion that Pc must 
be zero if < A;^ > is infinite, certainly in one dimension. In 2 dimensions the 
situation is complex, and will take more work to fully sort out. 
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